Particle creation due to tachyonic instability in relativistic stars 
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Dense enough compact objects were recently shown to lead to an exponentially fast increase of the 
vacuum energy density for some free scalar fields properly coupled to the spacetime curvature as a 
consequence of a tachyonic-like instability. Once the effect is triggered, the star energy density would 
be overwhelmed by the vacuum energy density in a few milliseconds. This demands that eventually 
geometry and field evolve to a new configuration to bring the vacuum back to a stationary regime. 
Here, we show that the vacuum fluctuations built up during the unstable epoch lead to particle 
creation in the final stationary state when the tachyonic instability ceases. The amount of created 
particles depends mostly on the duration of the unstable epoch and final stationary configuration, 
which are open issues at this point. We emphasize that the particle creation coming from the 
tachyonic instability will occur even in the adiabatic limit, where the spacetime geometry changes 
arbitrarily slowly, and therefore is quite distinct from the usual particle creation due to the change 
in the background geometry. 

PACS numbers: 04.40.Dg, 04.62. +v, 95.30.Sf 
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I. INTRODUCTION 

It was recently shown that relativistic stars may be- 
come unstable due to quantum- field effects P, Q|. The 
so-called vacuum awakening effect occurs for a free scalar 
field <f> properly coupled to the spacetime curvature 
This effect is characterized by an exponential point- 
dependent increase and decrease of the vacuum expecta- 
tion value of the stress-energy- momentum tensor (T a b) . 
This is caused by a tachyonic-like instability, which in- 
duces an exponential growth of (<I> 2 ). Once the effect is 
triggered and the scalar field exits its initial quiescent 
regime, few milliseconds would be enough for the vac- 
uum to overwhelm the energy density of the compact 
object. The star destiny is presently uncertain because 
it depends on how scalar field and spacetime geometry 
evolve in the unstable phase to reach a final stable config- 
uration. As recently argued in Ref. [3[ , in some cases the 
appearance of a proper scalar field could restabilize the 
star, a phenomenon usually called spontaneous scalariza- 
tion Q. This would typically change the star gravita- 
tional mass by a few percent. However, depending on 
how the star enters the unstable phase, it seems possible 
that the scalar field does not react fast enough, leading 
to some dramatic implosion/explosion event. Whatever 
turns out to be the final configuration, being the star 
somehow rebalanced or destroyed, the unstable phase 
must be detained and the vacuum must evolve to some 
new stationary regime. This observation alone allows 
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us to extract quite important information about the fi- 
nal state of the scalar field as the vacuum "falls asleep" 
again. In particular, we show that the vacuum fluctua- 
tion built up during the unstable epoch leads to particle 
creation in the final stationary state. The amount of cre- 
ated particles will depend mostly on the duration of the 
unstable epoch and final stationary configuration. 

The paper is organized as follows. In Sec. UTJ we dis- 
cuss the quantization procedure for a free scalar field in a 
curved spacetime in the presence of tachyonic-like modes. 
In Sec. IIII( we apply the previous-section results to re- 
view the vacuum awakening effect in relativistic stars. 
In Sec. IIV1 we probe the unstable phase using Unruh- 
DeWitt detectors. In this period, no natural particle 
content can be ascribed to the scalar field and the use of 
detectors is particularly useful to investigate the behav- 
ior of the vacuum fluctuation. We show in this section 
that even assuming a static spacetime in the unstable 
phase when the vacuum is "awake" , particle detectors 
following orbits of the timelike isometry will copiously 
excite. This is possible according to co-static observers 
because each detector excitation is accompanied by a cor- 
responding decrease in the energy stored in the field due 
to the excitation of a nonstationary (tachyonic) mode 
which contributes with negative energy. Then, in Sec.fVl 
we show that at least part of the quantum fluctuations 
built up in the awaken phase eventually draws to particle 
creation as the unstable period ends and the vacuum falls 
dormant again. We emphasize that the particle creation 
occurs even assuming that the spacetime change is arbi- 
trarily slow and, thus, differs, e.g., from the well-known 
phenomenon of particle creation in evolving universes in- 
duced by the change of the background geometry 
A toy model is also offered to illustrate in a concrete sce- 
nario the conclusions above. We close the paper with our 



final remarks in Sec. I VII 

II. PRELIMINARIES: FREE SCALAR FIELD 
QUANTIZATION WITH TACHYONIC-LIKE 
MODES 

A. Standard field quantization in globally 
hyperbolic spacetimes 

In this section, we review the standard quantization 
of free scalar fields in curved spaces [1, HJ giving par- 
ticular attention to the case of static spacetimes with 
tachyonic-likc modes. Let us begin by considering a glob- 
ally hyperbolic spacetime (Ai,g a b) foliated with Cauchy 
surfaces St labeled with a parameter t. Now, let us cover 
the S t surfaces with x l (i = 1, 2, 3) coordinates satisfying 
n a V a a; 1 = 0, where n a is the future-directed unit vector 
field normal to Et. In terms of coordinates x = (t,x l ), 
we can write the spacetime line element as 

ds 2 = N 2 {~dt 2 + hij(x) cfeW), (1) 

where N = N(x) > is the lapse function and ^gij — 
N 2 hij is the three-dimensional spatial metric induced on 
each Cauchy surface S t . 

We define the dynamics of a real scalar field $ with 
mass to in (Ai,g a b) through the action 

S=~ I d 4 2V^(V a $V a $ + 777 2 $ 2 +ei?$ 2 ), (2) 

where g = det (g a b) and £ G R determines the non- 
minimal coupling between the scalar field and the scalar 
curvature R. This leads to the following Klein-Gordon 
equation: 

-V a V a $ + m 2 $ + £i?$ = 0. (3) 

Next, we define the Klein-Gordon inner product between 
any two solutions u and v of Eq. Q as 

(u,v) KG = i I dT,n a [u*V a v-vV a u*}, (4) 

where dT, is the proper- volume element on St, and we 
recall that Eq. (Q} does not depend on the choice of E t . 
The conjugate-momentum density Tl(x) is defined as 

n = 6S/6& = „ Q v a $, (5) 

where " ' " = d t and ^g = det(( 3 ) gy) . The canonical 
quantization procedure consists of promoting field and 
momentum density to operators $ and II, respectively, 
satisfying canonical commutation relations: 

[*(t,x),$(f,xO] E( - [n(t,x),n(t,x')]E t =o, (6) 

[$(t,x),n(i,x')]s t - ^ 3 (x,x'), (7) 
where x = (x , x 2 , x 3 ). 



In order to realize a representation of these commuta- 
tion relations, consider positive- and negative-norm so- 
lutions of Eq. (|3]), u'a" 1 and u« "* = (u^)*, respectively, 
which together form a complete set of normal modes sat- 
isfying 

(ui+),4 +) ) KG - -(ui-\u { - ) ) KG = S(a,p), (8) 
(ui + \u ( ^) KG = 0. (9) 

Here, S(a, /3) is the delta function associated with the 
quantum numbers formally described by a, /3. Then, we 
construct the field operator using } as 

$> = J dn{a)[a a u^ +aW~\ (10) 

where /i is a measure defined on the set of quantum num- 
bers and in order to satisfy Eqs. ^ and ([7]), the annihila- 
tion a a and creation d\ operators must satisfy the usual 
commutation relations [a Q ,aJj] = S(a,f3), [a Q ,a,3] = 0. 
The vacuum state |0) of this representation is defined by 
requiring a CT |0) =0 for all a. 



B. Quantum fields in static spacetimes 

Now, we restrict our analysis to static spacetimes in 
which case the line element ([I} is cast as 

ds 2 = N 2 (-dt 2 + %(x) dx l dx j ), (11) 

where N = -/V(x) > 0. Under this condition, we write 
the field equation © in the form 

— = [-A + Ks(* (12) 

where $ = N$, A is the Laplace operator associated 
with hij, and 

^(x) = N^AN + N 2 (m 2 + £i?) 

= (1 - G^N^AN + N 2 m 2 + £K (13) 

is the effective potential with K = K(x) being the scalar 
curvature associated with hij . 

The existence of a timelike Killing field n b — (dt) b 
associated with spacetime (fTTj) suggests that we look for 

solutions of Eq. ([12")) in the form cx F a (x) exp(— iuj a t) 
corresponding to solutions oc N~ 1 F a (x) exp(— iuj a t) 
for Eq. ©. In this case, F a (x.) will satisfy 

[-A + V; ff (x)] i^(x) = A CTJ F CT (x), \ a — to 2 , (14) 

with proper boundary conditions. At this point, the only 
restriction on A CT is the one imposed by Hermiticity of the 
operator — A + V c ff(x), which demands X a G K. 
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Let us consider first solutions of Eq. (fT4|) with posi- 
tive eigenvalues: A CT = w 2 a > 0. Then, the correspond- 
ing positive-norm solutions satisfying Eq. (j3|) will be the 
usual oscillatory modes: 



,(+) 



-F a (x), zu a > 0, 



where we demand 

d 3 xVhF a {x)*F (x.) = 5{a,P) 



(15) 



(16) 



in order to guarantee that modes (|15j) are properly nor- 
malized according to Eqs. © and ©. 

Now, we note that in some cases Eq. (TT4)) also allows 
for solutions with negative eigenvalues: X a = — fi^ < 
0. These solutions are associated with solutions of 
Eq. (|12[) with exponentially increasing and decreasing 
exp(±f2 CT i) time dependence. Under such circumstances, 

{vi + \v& "*} must be supplemented by an extra set of 
modes {wi + \wi '} in order to generate a basis for the 
solution space of Eq. ([3]). Normalized positive- norm 



modes 



,(+) 



can be found and read 



n fi(r t— i 



w {+) = e iP„ 



— n a t+ia a 



yjAVt a sin(2a CT ) N(x) 



-f ct (x), n a > o, 



(17) 

where a a € ]0, j] and, for the sake of convenience, we 
did not vanish the arbitrary global phase f3 a yet. By 
choosing, e.g., a a = fi a = n/4, Eq. (fTT)) would assume 
the simple form 



(« 



2 v / fWV(x) 



F CT (x), n a > (18) 



but we shall adopt here the same choice as in Ref. [l|, 
where (3 a — and a a = 7r/12: 



,(+) 



V2n^v(x) 



F tr (x), > 0, 
(19) 



in order to make look "as similar as possible" to 

vi + \ Because and wi~^ grow exponentially in time, 
we borrow from cosmology the "tachyonic" term (see, 
e.g., Ref. [10]) and refer to these modes accordingly. 
(It should be noted, however, that in the cosmological 
context the scalar field is self-interacting, as ruled by 
some interacting potential, in contrast to our present case 
where it is free.) 

As a result, the field operator $(a;) can be constructed 
using VcT^} and {u4 + \w^r'} as 



(20) 



where [6 a , St] = £(a, /3), [c a , ct] = <5(a, /3) (with the other 
commutators vanishing), and we have used the summa- 
tion symbol in the right-hand side of Eq. (|2U|) because the 
tachyonic modes will be labeled later with quantum num- 
bers assuming discrete values. We recall that the vacuum 
state |0) satisfies b a \0) = c CT |0) = for all a. We note 
that in contrast to the v& and modes, and 
Wa are not frequency eigenstates of id/dt. As a result, 
the vacuum |0) and the other Fock states do not have in 
general any natural particle-content interpretation (see 
Refs. [11] and [l2[ for a field-theoretic discussion on the 
Fock space in the presence of tachyonic modes). This fea- 
ture will lead us to use Unruh-DeWitt detectors to probe 
vacuum fluctuations of the scalar field in Sec. IIVI 

Nevertheless, important pieces of information are di- 
rectly provided through the (formal) stress-energy tensor 
operator: 

f ab = (1 - 2£)V a $V b $ + £R ab $ 2 - 2£$V a V b <l 
+ (2£-l/2) 5ab [V c <IV c l> + (m 2 +eR)<l 2 ] (21) 
and the corresponding Hamiltonian: 



H 



dY> a Kb T 



ah 



dYi p, 



(22) 



where e?£„ 



dSn a , K b = (d t )b, 

P = n a K b T ab 



(23) 

is the energy-density operator in S t associated with the 
timelike isometry, and we recall that Eq. (|2"21 does not 
depend on the St choice because \7 a (nbT ab ) = 0. Thus, 
the total energy is conserved. By using Eq. (f2"U)) in the 
Hamiltonian operator (|22|) . we obtain 



H 



dY,N~ lr f< 



St 



= \J dix{a)(bth+b„bt)za a +Y,ticr, (24) 

where 

ii a = -[V3/2{cU* + c a ci) + c a c a + £gt]Or. (25) 
In contrast to the first term in the right-hand side of 
Eq. (|2"4"]l . associated with the oscillatory v& modes, 
which always gives a positive-definite contribution to the 
energy expectation value for every state choice, the sec- 
ond term, associated with the tachyonic modes w„ , 
gives a negative contribution to the energy expectation 
value for some states. This can be easily seen by rewrit- 
ing Eq. ([25]) as 

H a = (l-V3/2)pl-(l + V3/2)nlq 2 a , (26) 

where we have defined the position- and momentum-like 
operators 

q a = -JL^ + gt), (27) 



V2TT 



(28) 
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respectively, satisfying [q G , p a ] — i I with / being the 
identity operator. Equation (|26j) is formally identical to 
the Hamiltonian of a non-relativistic particle in a har- 
monic potential turned upside down [131 ] ■ It is clear, then, 
that the "potential" term gives a negative contribution 
to the energy expectation value. In particular, for states 
| <3>) satisfying 



|*>=3|¥>, SeN, 



(29) 



which include the vacuum state, it is easy to see that 
and, thus, 



(*|c ff c CT + ctct|*) 



(tf |7k|tf) = -v / 3(l/2 + 5) n a < 0. 



(30) 



Hence, for these states the negative contribution coming 
from the "potential" in Eq. (|2"o]l dominates over the cor- 
responding positive one coming from the "kinetic" term. 
(We shall return to this point when we discuss the exci- 
tation of Unruh-DeWitt detectors in Sec. IIV1 ) The fact 
that the right-hand side of Eq. (|3"0"]) may be arbitrarily 
negative for sufficiently large S reflects the fact that T-L a 
is unbounded from below. 



III. AWAKING THE VACUUM IN 
RELATIVISTIC STARS DUE TO TACHYONIC 
INSTABILITY 

Now, we shall see how tachyonic modes can appear in 
neutron-like stars and discuss their consequences 0, ■ 
Let us assume the case in which (A) classical matter ini- 
tially scattered throughout space with very low density 
eventually collapses to form (B) a static and stable star 
according to general relativity. Spacetimes associated 
with situations A-B are well described by the line ele- 
ments [see Eq. ([TT])] 



ds 2 = 



(A) 

iV ( 2 B) (~dt 2 + dx l dxi ) (B) 



-dt 2 

t2 



(31) 



We note that for the time being we will restrict our in- 
vestigation to the static regions A and B of the space- 
time. Comments about how our present analysis can be 
completed as one takes into account the time evolution 
between the static eras are made along the text. For 
the sake of obtaining explicit results, we make an ex- 
tra simplification in this section and consider spherically 
symmetric stars in which case Eq. (|3"Tj) is replaced by 



ds 2 



-dt 2 + dx 2 (A) 
-f(dt 2 - d X 2 ) +r 2 (d6 2 + sin 2 d dip 2 ) (B) ' 

(32) 

where f = fix) > ® an( l r = r (x) > satisfy fix) —> 1 
and r(x)/x —> 1 for x ~~ ^ °°! an d drjd\ > so that no 
trapped light-like surface is present. 

We construct the field operator similarly as in Eq. (fTDl : 



(33) 



where we choose here u^f* such that they assume the 
usual flat-space stationary form in the asymptotic past 
(region A): 

( = ] {Wn 3 ^)- 1 / 2 exp[TiKi - k • x)] (34) 



with k e R 3 and Wk = Vk 2 + to 2 . This choice is moti- 
vated by the fact that we shall assume hereafter the scalar 
field to be in the no-particle state |0)i n as described by 
static observers in the asymptotic past: dk|0)i n = 0. 

Now, let us represent $ in terms of positive- and 
negative-norm modes in region B, when the star has set- 
tled down, as [see Eq. (|20|) ] 



(35) 



where [3] 

,(+) (B) mt F^iix) 



v , = e 



(36) 



,(_+) @ / nt_i W /i2 



■ ^L (x) Y ltl (e,i>), n>o. (37) 



'2fir( X ) 

Here, 5^(0,0) are the usual spherical harmonics (I = 
0, 1,2, . . . and n = -1,-1 + 1,.. .,1), F m iix) and Fm(x) 



satisfy 



and 



J^ F ™ + V^Fm = -n 2 Fm, 



respectively, and 



V. 



(0 



off 



/ [m 2 +CR 



1(1 + 1) 



1 d 2 r 



rdx 2 



(38) 



(39) 



(40) 



is the effective potential. For perfect-fluid stars, the ef- 
fective potential (|40"]) can be cast as 



vS = f 



2 , ^ + 1) 

TO "I n 



1\ 8nG, 



(41) 

where p = p(x) is the mass-energy density of the stellar 
fluid and 



pix) 



3M( X ) 
47r[r(x)] E 



(42) 
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is the corresponding average density up to the radius co- 
ordinate r{x)i which encompasses a mass M(x)- We re- 
mind that according to general relativity R = 8nG(p — 
3p), where p is the hydrostatic pressure which bears the 
star up against its weight. 

As discussed in Sec. Ill Bt the appearance of tachyonic 
modes in the present context will depend on the existence 
of nontrivial solutions for Eq. (13T))) . They are expected to 
exist for negative enough effective potentials satisfying 



(43) 



where r = R s is the star radius. Because the centrifugal 
barrier in Eq. (|4"Tj) is positive, we look for tachyonic so- 
lutions of Eq. (|39|) with I — which are the most likely 
ones to exist (if any). By taking / ~ 1 and assuming 
p ~ p, we obtain from Eq. (|4"Tj) that 

V$ ~ m 2 + (£ - 1/6)71. 

Clearly, only the second term in the right-hand side can 
be negative. Then, Eq. (|43p implies that a necessary 
condition for the existence of tachyonic modes with £ rs 1 
is 



10 15 g/c 



7 km 



and 



m 2 /(3.5 x lfT 11 eV) 2 
p/(10 15 g/cm 3 ) 



> 1 



« 1, 



(44) 



(45) 



where we have set p ~ p. Equations (1441) and (|45[) show 
that the appearance of tachyonic modes for small £ val- 
ues in the spacetime of typical neutron stars requires the 
scalar field to be light: m < 1CP 11 eV. 

Although light scalars are widely considered in astro- 
physics and cosmology, there is the issue about how much 
extra mass they could acquire from Planck-scale radia- 
tive corrections. For axions, e.g., a general expression for 
the mass shift can be cast as Sm 2 ~ K a fa +2 /Mp, where 
K a is some unknown coupling constant, f a ~ 10 12 GeV 
is the energy scale of the Peccei-Quinn symmetry break- 
ing, Mp ~ 10 19 GeV is the Planck energy, and n is a 
mo del- dependent positive integer (associated with the di- 
mension of the symmetry-breaking operators appearing 
in the effective Lagrangian) We see, then, that Sm a 
can easily exceed, e.g., 10 -5 eV (ruling out axions as a 
dark matter candidate) unless K a and n turn out to be 
small and large enough, respectively (see Ref. [l6| and 
references therein). Fortunately, explicit models showing 
how scalar fields can be protected from acquiring large 
mass due to quantum gravity effects have already been 
worked out (see, e.g., Ref. |17|)- I n our context, assum- 
ing the electroweak symmetry breaking of the standard 
model which has an energy scale of Aesm ~ 100 GeV, the 
corresponding mass shift would be 6m 2 ~ KA^/M™, 
where again K and n are unknown. Here, we pragmat- 
ically assume that Planck-scale effects will not shift the 



mass of our scalar field beyond 10 _12 eV. At this point, 
it is difficult to say how strong this assumption is because 
of our lack of understanding of the Planck-scale physics. 
Still, this is much less demanding than what is usually 
required for quintessence fields where the mass shift can- 
not typically exceed the mass scale defined by the Hubble 
constant 10~ 33 eV [llj]. A detailed analysis of this issue 
would be welcome but goes far beyond the scope of this 
paper. For computational purposes, we take our scalar 
field to be massless. 

By assuming stars with uniform and parabolic den- 
sity profiles and suitable £ values (typically £ > 1/6 and 
£ < —2), it was shown in Ref. @ that tachyonic modes 
do appear for M/R s ratios compatible with neutron-like 
stars. 

Now, let us proceed by recalling that the positive-norm 
in- modes it^ , which in region A look like as exhibited 
in Eq. (|34[) , will emerge, in general, as a combination 

of positive- and negative-norm modes {v^^v^^} and 



..(+) w (-) 



t } in region B [see Eqs. (l36l) and (|37|) ]. Hence, 



not only the in-vacuum will not coincide in general with 
the out-vacuum but also at least some of the in-modes 
will certainly go through a phase of exponential growth 
provided, of course, the existence of tachyonic modes 
tUfj^. This leads to what was denominated vacuum awak- 
ening effect in relativistic stars, i.e., an exponential am- 
plification of the vacuum fluctuations [H, Hj] . I n order to 
see this, we use Eq. (1331) to calculate 



|l> 2 |0)i 



(B) e 

~ K 



2S2 ( 



8irti\ r( X ) 



[l + 0(e-*% (46) 



where Fq (x) denotes the solution of Eq. (|3U)) with the 
most negative eigenvalue, — fi 2 (taking I = 0, which is the 
most favorable case), e = const > 0, and k is a constant 
of order unity whose value depends on (i) projections 

of modes on and (ii) the quantum state, as- 

sumed here to be the in-vacuum |0)i n . It is worthwhile 
to emphasize that ultraviolet divergences, which should 
be renormalized to obtain ($ 2 ), are associated with the 
ki^cxd sector of the oscillatory modes [see Eq. (|35"]) ] and 
does not concern the tachyonic modes (fi 2 < fr < oo), 
which are the ones giving the dominant contribution in 
Eq. (|46[) (because of the exp(2f2i) term). Accordingly, 
the expectation value of the vacuum energy density (|23[) , 
namely, 



»(0|p|0)i 



i(0|T a6 |0)i 



experiences an exponential amplification: 
- D 2fit [i_4£ d ( 



(b) f2e 



2 d ( F no 



1671-7/ I 2r 2 d X \ dx V Mr 



j d (F^df 
H 2 r 2 d X \ f d X . 



■[l + 0(e- et )]. (47) 



The time scale which rules how fast the vacuum energy 
density increases is given by f2 _1 ~ | -1 / 2 ~ R s [see 
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Eq. (|43)) ]. By using this, we rewrite Eq. (|47j) as 



1 (0|/5|0) i 



(B) 



(B) 



nh(f)e 2nt /R 



h(r) exp 



t/(10" 5 s) 
i? s /(10km) 



10" 62 g/cm J 
i? 4 /(10km) 4 ' 

(48) 



where h(f) is a dimensionless function of f = r/i? s (which 
vanishes asymptotically and is of order unity for f ~ 1). 
We see from Eq. (f48|) that once the effect is triggered by 
a neutron-like star with R s 10 km, few milliseconds 
would be enough for the vacuum energy density to be- 
come dominant over the star classical mass-energy den- 
sity (which can be as high as 10 14 — 10 17 g/cm 3 ). We must 
emphasize that at some point the spacetime must back- 
react against the growth of the vacuum energy density, 
affecting the field and ceasing the instability by taming 
the tachyonic modes. Eventually, field and spacetime 
must reach a new stable configuration. As argued in 
Ref . [|| , one possibility would be that for some values of £ 
the spontaneous scalarization mechanism [3] could resta- 
bilize the star. In our context, this would correspond 
to a symmetry breaking which would lead ($), which is 
null as calculated in the |0)i n vacuum state, to acquire 
a nonzero large value compatible with the exponentially 
amplified ($ 2 ) [see Eq. (|4"o) ]. Whether the star will end 
up destroyed or somehow rebalanced is unknown at this 
point. 

We close this section explaining how the vacuum en- 
ergy density amplification is consistent with energy con- 
servation discussed below Eq. ([25)1 . For this purpose, let 
us note that V a (K&T ab ) = can be rewritten as 



d t p- 



(49) 



fl 



where j = — ^/JT'q is the energy-current density. The 
corresponding vacuum expectation value can be calcu- 
lated and reads 



•|J |0)i 



(B) 



-nV l (f)e 2Qt /R 



(50) 



where ^V z V l = g)- 1/2 d l {^/^ V 1 ) = 2Slh(f). 
Thus, the total energy is conserved because the gravi- 
tational field redistributes the vacuum energy density in 
such a way that an amplification of i n (0|/5|0); n somewhere 
with positive magnitude must be compensated elsewhere 
by a corresponding amplification with negative magni- 
tude. 



IV. PROBING THE AWOKEN PHASE USING 
DETECTORS 

Now, in order to probe the building up of the vac- 
uum energy density in region B, where the vacuum is 
awake by the presence of tachyonic modes, we will look 



directly at the response of Unruh-DeWitt detectors. We 
shall do so because, as discussed in Sec. Ill BL the Back- 
space states have no natural particle-content interpreta- 
tion [see discussion below Eq. (|20[)]. Here, we relax the 
spherical symmetry assumption of the previous section 
and consider regions A and B as described by the line 
elements (I31[) . Because we want to avoid any contribu- 
tions in the response coming from the motion of the ap- 
paratus, the detector (with proper time r) is made to lie 
static following an integral curve x = x(t) of the timclike 
isometry. 

We consider here a two-level Unruh-DeWitt detec- 
tor represented by a Hermitian operator too acting in 
a Hilbert space spanned by unexcited and excited energy 
eigenstates \E ) and \E) (E > E ), respectively. The 
detector is prepared to be initially unexcited. For our 
purposes, it is convenient to switch it on in the very be- 
ginning of region B, where we set r = 0. At the tree level, 
the excitation probability as a function of the proper-time 
interval T is given by 



Lity 

[7] 



P cxc = \(E\m \E )\ 2 T{AE), 
where AE = E — Eg and the response function is 

i-T pT 



(51) 



T{AE)=j dr f dr'e- iAE ^-^G+[x(r),x(T')] 
Jo Jo 



(52) 



with 



GiN(r), a ;(r')]^ i „(0|$[x(r)]<I.[a ; (r')]|0) i 



The two-point function is written throughout the space- 
time in terms of the in-modes as 



= / d 6 k 



4 +) (aO«L _ V), 



(53) 



where we recall that in the asymptotic past u^\x) take 
the simple form given in Eq. (|34[) . Next, we suitably 

decompose w k + ' in terms of {«S } and {u^} as 

(+) * (+) a (-) 

+ / ^(a)« k ^+)-^ k ^-)) ! (54) 



where, for the sake of simplicity, we have assumed that 
the scalar field is made unstable in region B by the exis- 
tence of a single tachyonic mode [see Eq. (|19p]: 

(Mt-i-n 111 I p -nt+i-rr/12\ 

.„(+>© (e '_+e M^) (x)j Q>Q 



(55) 

The tachyonic mode above will be labeled by f2 (with 
the other quantum numbers being omitted to simplify 
notation), while the oscillatory modes are cast here as 
[see Eq. §15$} 



,(+) <!> 



%/2^JV (B )(x) 



^i B) (x), tn CT >0. 



(56) 



G 



The Bogoliubov coefficients in Eq. (|54|) are 

a<rk = (4 +) ."i +) )KG, Ark = ,«£ +) )kG) 

ank = («k '^n )kGi /?ok = -(4 _) ' w !! +) )kg 

as calculated in any Cauchy surface. 

Now, we must proceed and evaluate the response func- 
tion (|52|) . For this purpose, it is convenient to use 
Eqs. ([5i ]) -(|56 ]l to calculate 



dTe- lAET u { +\x(T)] = / dfx(a) 



F 



(B) 



(B) 



-tp a+ ~ /3, 



F 



(B) 



<rk 



(B) 



,,(B) 



? (B)* 



* ±n_ 

*Ok 



(B) 



Ok 



(B) 



(57) 



where we have defined 



dr- 



^ e nr/JV CB )-i7r/12 + e -Or/Af (B) +i7r/12^ 



2fl 



e Wl2(l _ g-(fi/iV (B) ±iA.E)T\ 

n/AT (B) ± iAE 



(58) 



Ti = \ d A k 



d/J,(a) 



F 



(B) 



(B) 



-if) a+ - /3, 



(B)* 



ok 



AT, 



(B) 



The physical meaning of the response J-(AE) is more 
easily grasped in the case where the proper time interval 
T is "large", i.e., T > AiT 1 (in addition to T > rr 1 
whenever the tachyonic mode is present) and we will as- 
sume this hereafter up to the end of this section. In this 
case, Eqs. (|5"8"j) and (|59")l can be written as 



(64) 



9- 



-i(7r/12±A_E)T TO/AT (B ) 



and 



'27T 2 



-i(a„/Nm ) ±AE)T/2 



S(m a /N {B) ± AS), 



(65) 



respectively. For the sake of comparison, we shall discuss 
separately the situations where the tachyonic mode 
is present from the one where it is absent. 

In the case where the tachyonic mode Wq is absent, 

(B) 



ply 



and, thus, the detector response becomes sim- 



T(AE) = J- 2 | x=Xd . 



l 



dT e - l (^/JV (B )±AS)r 



\/2w a JO 

2e -i(™ a /N (B) ±AE)T/2 



sin [(w a /N {B) ± AE)T/2] 
vo a /N {B) ± AE 



(59) 



and x = x<j is the detector's spatial position. Then, we 
write the detector response function ((52)) with the help 
of Eq. (J57|) as 



F{AE) = [To + F,+ F 2 ] y 



(60) 



This is the usual result whose interpretation is straight- 
forward: assuming that the detector stays switched on for 
an arbitrarily long time T, we have from Eq. (|65[) that 
« and according to Eq. (|63j) the detector excites 
by the absorption of particles created due to the space- 
time transition from regions A to B (/3 a u ^ 0). (Note 
that when no restriction is posed on T, the detector ex- 
citation will also have a contribution coming from the 
process of switching it on and off (ipo-+ 7^ 0) [l9j|.) 

It is also interesting to note that in the absence of the 
tachyonic mode, the response will not grow faster than 
T. By recalling from Eq. (|65|) that ip a + R3 0, we write 
the response as [see Eq. (l63|) ] 



where 



d A k 



(B) 



rv* Fq 



- A 



F 



(B)* 



Ok' 



N, 



(B) 




F 



(B)* 



c+ - a 



(B) 



<rk 



AT, 



(B) 



(61) 



(62) 



« . 771(B)* 
J-(AS) » J d A k J M°)j^ 



"A<rkVC- 



< / d A k[ / dMWI^kp 



x / dM(o-)|^<7-| 



2 I^ B) I 2 



(66) 



where we have used above the Cauchy-Schwarz inequal- 
ity. Now, by using Eq. (|59[) [or directly Eq. ([65]) ] in 
conjunction with the identity limA-^oo sin 2 (oj A) /uj 2 A — 
tt5(lu), we obtain 



{■nT /m„) 5{m a /N(p) - AE). 



(67) 
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This is used in Eq. (|66j) to conclude indeed that the re- 
sponse will not grow faster than T: 



F{AE) < C 2 T 



(68) 



with 



C 2 



d 3 k[ / d/i(a)|/? CTk | : 



d[i(<j) 



\F t 



(B)|2 



"(B) 



8{uj a - N (B) AE) . (69) 



On the other hand, assuming that the tachyonic mode 
Wq is present, the detector response will be dominated 
by Eq. (jCTj) and, thus, F(AE) « J r o| x=Xti - Thus, we use 
Eq. dM]) in Eq. dH]) to obtain 



T(AE) « Z exp[2ra/Ar (B) (x d )] 



x / d 3 fc 
where 



i 7 , 



(B) 



£2 -in/12 



A', 



(B) 



(B) 



Qk 



£2 „W 12 



A', 



(B) 



(70) 

x=x d 



2n[AE 2 + nyNf B) { Xd )Y 

The exponential increase in the detector response reflects 
the growth of the vacuum fluctuations and will continue 
as long as the unstable phase is not forced to terminate. 
This is possible because each excitation of the detector is 
accompanied by a corresponding decrease of the energy 
stored in the field due to the excitation of a tachyonic 
mode Cfj|0)in with negative energy expectation value [see 
Eq. (f5U)) ] and corresponding discussion). The copious 
excitation of the detector realizes the fact that in the 
unstable phase the scalar field functions as an energy 
reservoir only limited by backreaction effects. 



V. FALLING ASLEEP OF THE VACUUM AND 
PARTICLE CREATION 

A. General discussion 

As already mentioned, at some point the unstable 
phase must cease, leading the system back to some sta- 
tionary configuration. This will be represented by the 
static region C in Eq (|TTj) . which completes the scenario 
presented by Eq. (pTTj) : 



-dt 2 



ds 2 



rfx 2 

+2 



(A) 



N (B) i- dt ' 2 + h^'dxHx*) (B) 



(71) 



N 2 c) (-dt 2 + h>^>dx l dx3) (C) 



Here, iV(j) = iV(j)(x) > 0, J e {B, C}, are smooth func- 
tions and h^p — h\p(x) — 1, 2, 3). We note that for 



the sake of simplicity we are using the same coordinate 
notation (t, x) for the three epochs. 

Because the field is assumed to be deprived of tachy- 
onic modes in region C, we expand $ as 



(72) 



where {v„ \ Va } are normal modes which in region C 
satisfy 



(+) (°) e " 



V2^7AT (c) (x) 



Ff)(x), m a >Q, (73) 



and analogously for ui ' . All symbols in Eq. (l73|) can be 
inferred from Eq. (|56"|) by replacing "B" by "C". More- 
over, because region C is also static, we may wonder what 
will be the particle content of the scalar field in this re- 
gion. The key point consists in realizing that the in- 
vacuum fluctuations which were exponentially amplified 
during the unstable phase cannot, in general, be accom- 
modated as mere fluctuations of the out-vacuum state 
(the one defined by rf CT |0)out = for all a and which rep- 
resents absence of particles according to static observers 
in region C). In conclusion, a burst of particles is ex- 
pected as the field exits the unstable phase B. 

Let us estimate the expectation number of created par- 
ticles in the simplified case where the spacetime is sym- 
metric by time reflection with respect to some Cauchy 
surface S fs in region B. Hence, we consider a particular 
case of Eq. (|TTj) . namely, 



-dt 2 



ds 2 = 



-dx 2 
-dt 1 ^ 



N 2 (-dt 2 + h\J ] dx*dx3) (B) , 
-dt 2 + dy? 



(A) 



(74) 



(C) 



although we emphasize that we have chosen regions A 
and C to be flat only for the sake of simplicity; the same 
reasoning presented below may be straightforwardly ap- 
plied to other static spacetimes. Now, we focus on 
the normal modes and i/£ with respect to which 
asymptotic observers in regions A and C define their no- 
particle states, respectively. In the past and future re- 
gions, they assume the following forms: 



,(±) (A) 

■ (C) 



(167r 3 w k )- 1/2 exp[=Fi(w k < - k • x)], (75) 



(167r a w k )- 1/2 exp[T«(w k i - k • x)]. (76) 



We are interested in an d ^ k evolved forward and 
backward to the beginning and end of region B, namely, 

(to,x) and (to + ?\ x )> respectively. Here, t = 
to = ts — T/2 determines the beginning of the unstable 
phase B and T represents its coordinate-time duration. 
Because of our assumption that the spacetime is sym- 
metric by time reflection with respect to Y* ts , we have 
(up to global phases) 



,(+) 
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In particular, for t = T/2: 

^(to + T, X ) { = } ui +) (t Q , X ). 

Then, we use Eq. (|54| to decompose in terms of i4 
and Wq^ , obtaining 

v -k (*0 + T : x ) = «?2k u; f2 +) (*o, x) - fck4~ } ( 4 o, x) 



(77) 



dM(<T)[a* k «^ +) (to,x) 
/3<rk^ _) (io,x)], 



where we assume again the existence of a single tachyonic 
mode for the sake of simplicity. 

In order to investigate particle creation in region C, 
we must, e.g., project it k (to + T, x) into i/£ '(to+T^x), 
where [see Eq. ([54)) ]: 



?4 +) (to+T,x) = ank«4 (*Q + T i x ) 

- /3nkWn~ Vo + T,x) 
d/i(a)K k 4+)(t + T,x) 

- /W^o+^x)]. (78) 

It is easy to see that 

i4±) (t + T, x) = exp(Tin7 ff T)^ ± ) (t„, x), (79) 
while we obtain from Eq. that 

w^Oo +T,x) = ±2isinh(^T=F«7r/6)w[i t) ^o,x) 

T2isinh(^r)w[ 2 =F) (to,x). (80) 

Then, by using Eqs. ([77]) and ([751) . we obtain for large 
enough f2T that 



\ u k> » u k )kg 



e Ckk' 



where the Klein-Gordon inner product was realized on 
the S( +t Cauchy surface and 

Ckk' = l [(af 2k e—/ 6 -/3ok)ao-k'-(«? 2k -/3 f2 ke M/6 )^_ k ,]. 

This leads to an expectation number of created particles 
with quantum numbers k' given by 

(Nv)~e™ T /d 3 fc|Ckk'| 2 , 



which grows exponentially as scaled by the product Q.T. 
In particular, even if the transitions from regions A to B 
and from regions B to C were made arbitrarily slow in or- 
der to minimize any particle creation due to background 
change (/3 CTk w 0), this would not alter the fact that a 
large amount of particles would be eventually created as 
the vacuum falls asleep (at least in the present scenario; 
see additional comments at the end of Sec. IV Bp . 

Next, we shall show that the burst of particles calcu- 
lated above does not rely on phase B being static; it will 
occur as long as the in-vacuum fluctuations get signifi- 
cantly amplified. 



B. A toy model 

In order to illustrate our general conclusion above, let 
us make an explicit calculation assuming a concrete sce- 
nario complying with the asymptotic static regions A 
and C considered in Eq. (|74|) but assuming some time 
evolution in the intermediate region B. This is in agree- 
ment with the idealized situation where initially spread 
out matter collapses to form a compact object and even- 
tually disperses back to infinity. Instead of calculating 
the particle production over the whole space, we shall 
restrict attention to the interior of a small cubical box 
with coordinate volume L 3 (oblivious to the matter form- 
ing the star), initially empty (of $ particles), placed in 
the very beginning at the spatial position where the star 
core will form. The convenience of introducing a small 
box is that we can cover its interior with approximately 
Cartesian spatial coordinates x (in which first-order spa- 
tial derivatives of the metric are negligible), writing the 
line element as 



ds 2 « a 2 (-dt 2 



dSi 2 ) 



(81) 



Here, a — a(t) > is introduced to reflect the background 
time evolution at the star's center (a = 1 in regions A 
and C) and we have omitted the second-order spatial 
dependence of the metric (which, nevertheless, contribute 
to the scalar-curvature term). The background evolution 
is chosen such that at some point the vacuum in the box is 
awaken by the presence of (six-fold degenerate) tachyonic 
modes [20]. After the unstable phase is terminated, we 
calculate the number of massless scalar particles which 
were created inside the box. 

Using Eq. (|8l]). we write Eq. ^ as 



,9£ 
~dt 



a z 



£m = o, 



(82) 



where V 2 = ^2-d 2 /dx^ is the usual Laplace operator. 
Assuming, for the sake of simplicity, periodic boundary 
conditions, we look for solutions of Eq. (1821) in the form 



</> k (t,x) 



Xk(t) 



^ik-x 



where k = 27rn/L with n e Z 3 (n 7^ 0) 
Eq. in Eq. (JH2J) we find that 



where 



dt 2 



Kff(t) 



" VeB(t) 



a 2 ZR 



Xk = k 2 Xk, 



-1j2 



d a/dt 



(83) 
By using 

(84) 
(85) 



Equation (|83p makes explicit another neat feature of in- 
troducing the small box: the boundary condition which 
it imposes locks the spatial dependence of the modes so 
that the time evolution can only mix modes with the 




FIG. 1: We plot —V e g as a function of t for Vb = 
1.6/(L/27r) 2 with the three smallest ojk = 27r||n||/L possi- 
ble values (see horizontal dashed lines). We also plot |Xk + '| 2 
for ||k|| = Lj£ in assuming Vo = 5L/2tt. Initially |x k +) | 2 
equals L/4n, then grows exponentially in the unstable region, 
where — V c g — (cj£" n ) 2 > 0, and eventually oscillates around 
(|/9k| 2 + l/2)(L/27r). The large amplitude which characterizes 
lxL + ^| 2 a t the end of the unstable phase reflects the fact that 
the in-vacuum fluctuations do not evolve into mere fluctua- 
tions of the out-vacuum state. 



same x dependence. This property will be used later to 
simplify the Bogoliubov-coefficient calculation. 

Now, we assume that energy density and pressure of 
ordinary matter at the center of the star drives R in 
Eq. (|85l) to induce the following simple form for the ef- 
fective potential: 



V eS (t) 



for t < and t > r/o 

4Vb - 1) for 0<t<77o 



, (86) 



where 770 , Vb = const > 0. We see that V e e(t) has a 
parabolic form in the region < t < t]q and reaches its 
minimum, — Vb, at t — 770/2 (see Fig. []}. 

Convenient sets of in-modes {ft } for i < and out- 
modes {V^ ± ' ) } for t > rjo are exhibited below: 



(±) *<o e 



^fi(uj^t — k-x) 



v. 



(±) t>»?o e 



=pz(wki — k-x) 



(87) 



where Wk = ||k||. The general expression of Ul ±] which 
complies with Eq. (|83|) and fits with its form (1571) in re- 
gion A is 



x) 



x k ±} W 



a(t)VL 3 



with < 0) 



/y / 2ui c . From the spatial de- 



pendence of the modes, we readily see that [recall the 
discussion below Eq. (|85|) ] 



ui + \t,Z) =a k V k (+) (*,x)+/?_ k vV(i,x), (89) 



where the Bogoliubov coefficients between the bases 
{Ul ±} } and {V^} are 



«kk' = ak"5kk', /?kk' = /?k"5k-k' 



For x k +) in region C, Eqs. (gTJ, ([55]), and ([89]) imply 



xL +) W 



t>r/ 



a iuj k t 



from where o?k and /3k = /3-k can be easily obtained in 
terms of x k + ^ an d Xk~^ = <^Xk evolved into region C: 



«k 



(-kx k +) + *xL +) 



Ac = 



(,, 



■ •(+) 

*x k 



t>Vo 



t>T)0 



(It can be easily verified using Eqs. (|84p and ((55)) that the 
expressions for a k and /3k above do not depend on the 
value of t > t/q.) Therefore, assuming that the field is ini- 
tially in the no-particle state |0)i n with respect to asymp- 
totic past observers as defined by the in-modes , the 



expectation value of created particles in region C [7ll21l 



|^out|0) in = 1/3 kk' 



k.k' 



is given by 

in<0|JV out |0)i 



k 

E 



K- / (+)|2 , , |,,(+)|2 
IX k I . w k|X k I 



2w k 



(90) 

t>T) 



where N out = £ k ^ utf d£ ut with d° ut and d° utt being the 
annihilation and creation operators defined with respect 
to the out-modes , respectively. This provides an 
expression for obtaining the expectation number |/3k| 2 
of created particles with quantum numbers k once the 
oscillatory in-mode Xk~^ ^ s (numerically) evolved until 
region C. 

The assumption to illustrate the appearance of tachy- 
onic modes consists in choosing a star which becomes 
dense enough and a coupling £ such that —V e ff — 
(a;™) 2 > for the least energetic (six-fold degenerate) 
modes allowed in the box, w£" n = 2tt/L, for some time 

interval. As a result, the corresponding xl solutions 
satisfying Eq. (|84|) are verified to exponentially grow for 
some time rather than to oscillate, triggering the vacuum 
awakening effect (see Fig.QJ. In Fig. [21 we plot |/3k| 2 as a 
function of 770, which scales with the time interval during 
which the vacuum stays awakened. 

Clearly, the final state is dominated by modes with 
u;£ lln = 2tt/L, which have experienced a phase of ex- 
ponential growth. The intensity of the particle burst 
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FIG. 2: The expectation value of created particles |/3k| 2 with 
quantum numbers k = 2irn/L is exhibited as a function of 770, 
where we have assumed that for some time interval the star 
becomes dense enough such that — V c ff — (^J? 111 ) 2 > 0, while 
for Wk > o;™ ln we always have — Ves — cu'i < (see Fig. [TJ. 
Here, we have chosen Vo(L/2iv) 2 = 1.6. 



sponding stress-energy tensor are 
1 



S = ~l d 4 x^(d a $d a <Z>* + V cB $$*) (92) 



and 



T ab = d {a <S>d b) <S>* - 



(93) 



respectively. Next, it can be shown from Eqs. ([91 
and (E! that 



da (T a b (d t ) b ) 



2 dt 11 ' 



which can be rewritten as 



dt J 2 dt 1 1 ' 



(94) 



(95) 



where 



1 / 9$ <9$* 



2 V dt dt 



and 



is strongly influenced by how long the vacuum remains 
awake. The inset of Fig. [2] focuses on modes with uj^ > 
2tt I L, which are not exponentially enhanced, and stresses 
the usually modest particle creation observed in time- 
varying spacetimes with asymptotic fiat regions @, [H[ . 
We note that in the adiabatic limit, where the back- 
ground geometry changes arbitrarily slowly (770 — > 00), 
particle creation for modes with Wk > 27r/L goes to zero 
as expected, in contrast to the ones for cjJJ 1111 = 2tt/L, 
which diverges. For a 10 m side box, an awakening time 
interval corresponding to 770 ~ 1CP 6 s would eventually 
lead to a massive creation of particles, with energy 2tt / L, 
engendering densities of 10 14 g/cm 3 , which is the typical 
density for some compact stars. If we relax our small box 
assumption and take L ~ 10 km, the same density would 
be reached for 770 ~ 10~ 3 s. Interestingly enough, this 
corresponds to the time interval for the vacuum energy 
density to take control over the evolution of the compact 
star once the vacuum awakening effect is triggered [see 
discussion below Eq. (|48|) ]. 

We stress that the conclusions above are derived by 
assuming the effective potential ([5rl)l . which is symmetric 
by time reflection, and can change depending on the star 
evolution. In order to show this, let us discuss the ener- 
getics of particle creation in the context of our toy model. 
For this purpose, it is useful to make the transformation 
<F -> $ = a$ to write Eq. (J52J) for $ as 



(91) 



[d 2 /dt 2 - V 2 + V eff (i)] 1 = 0. 



Thus, we have translated the problem into the simpler 
one of a scalar field $ in a flat spacetime (R 4 , 77^) subject 
to an external time-dependent potential V c g(t). 

The action which gives rise to Eq. (lint and its corre- 



1 d§*~ <9$ ~ 
j = - I ^t-V* + — V$* 



2 \ dt 



dt 



We see from Eq. (|95|) that the energy stored in the 
scalar field is not locally conserved whenever dV e g/dt 7^ 
0. The extra energy pumped into or out of the field is 
accounted by the "external agent" responsible to change 
V e g. Moreover, notice from Eq. (T9"51) that even if V e g(t) is 
symmetric under time reflection, the decrease in the field 
energy when it enters the unstable phase (dV e s/dt < 
and small vacuum fluctuations ($ 2 )) is more than com- 
pensated by the increase in the field energy when it exits 
the unstable phase (dV c s I dt > and large vacuum fluc- 
tuations ($ 2 ))- In fact, the latter can be overwhelmingly 
larger than the former, with the net extra energy being 
responsible for the particle burst. This analysis implies 
that in a physical situation, the final verdict concerning 
the amount of created particles will depend on a more 
detailed understanding on the spacetime evolution in the 
unstable phase, which would inform us about how long 
the vacuum would stay awake, and on the final classi- 
cal configuration reached by the gravitational and scalar 
fields, which would tell us how much energy would turn 
out available for particle creation. 

A closing remark for this section is in order. In our 
calculations the expectation value of the field remains 
zero throughout the background evolution. However, this 
field configuration ($) = is obviously unstable during 
the intermediate phase when the vacuum is awake. Thus, 
one may speculate whether during the transition to the 
intermediate phase the classical field profile (<&) could 
continuously adjust itself to nonzero values ("continu- 
ous spontaneous scalarization" ) in order to stabilize the 
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system, in which case tachyonic-like modes would never 
really be present. Unfortunately, a definite verdict to 
this question is beyond the scope of semiclassical gravity 
since it involves the subtleties of decoherence of a free 
field, initially in a state which is symmetric by the ex- 
change $ f> to a symmetry-broken phase in which 
(<f>) 7^0. Notwithstanding, a reasonable conjecture seems 
to be that coherence can be sustained for as long as the 
background geometry (which can be regarded as the sole 
classical "apparatus" with which the field interacts) is 
oblivious to <£>. In other words, it seems quite possible 
that (<&) = until backreaction becomes important. But 
when that happens, the fluctuations (<& 2 ) will already 
be amplified to the point where they cannot be accom- 
modated as mere vacuum fluctuations. A burst of par- 
ticles should follow, regardless whether ($) remains null 
or spontaneous scalarization takes place. Another impor- 
tant point is that spontaneous scalarization takes place 
only for negative values of £ [3(. Therefore, for £ > the 
whole scenario of a gradually changing ($) ensuring the 
stability of the system seems even more unlikely. 

VI. FINAL REMARKS 

After a review of the vacuum awakening effect in rela- 
tivistic stars, we have probed the exponential increase of 
the quantum field fluctuations using Unruh-DeWitt de- 
tectors. The fast increase of these fluctuations may lead 
eventually to an important burst of free-field particles af- 



ter the vacuum is forced to fall asleep again. This burst 
of particles would draw a significant amount of energy 
from the initial system. The amount of created particles 
will depend on the duration of the unstable epoch and 
on the final gravitational and scalar field configuration, 
which are open issues at this point. A possible signal fa- 
voring the vacuum awakening effect for a free field would 
be the unveiling of astrophysical events outpouring less 
amounts of visible energy than would be expected. This 
may also provide an efficient way of converting energy 
initially stored in the form of ordinary matter (forming 
the star) into a "dark" component which couples only to 
gravity. 

Simultaneous to the completion of this article, it was 
posted a classical analysis showing that during the scalar- 
ization process a strong emission of scalar radiation 
should occur [22| . which is in line with the conclusions 
presented here, especially with the discussion at the end 
of the previous section. 
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